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Space launch vehicles incorporate upper-level wind assessments to determine wind 
effects on the vehicle and for a commit to launch decision. These assessments make 
use of wind profiles measured hours prior to launch and may not represent the 
actual wind the vehicle will fly through. Uncertainty in the winds over the time 
period between the assessment and launch introduces uncertainty in assessment of 
vehicle controllability and structural integrity that must he accounted for to 
ensure launch safety. Temporal wind pairs are used in engineering development of 
allowances to mitigate uncertainty. Five sets of temporal wind pairs at various times 
(0.75, 1.5, 2, 3 and 4-hrs) at the United States Air Force Eastern Range and Western 
Range, as well as the National Aeronautics and Space Administration’s Wallops 
Flight Facility are developed for use in upper-level wind assessments on vehicle 
performance. Historical databases are compiled from balloon-based and vertically 
pointing Doppler radar wind profiler systems. Various automated and manual 
quality control procedures are used to remove unacceptable profiles. Statistical 
analyses on the resultant wind pairs from each site are performed to determine if 
the observed extreme wind changes in the sample pairs are representative of 
extreme temporal wind change. Wind change samples in the Eastern Range and 
Western Range databases characterize extreme wind change. However, the small 
sample sizes in the Wallops Flight Facility databases yield low confidence that the 
sample population characterizes extreme wind change that could occur. 


I, Introduction 

Space launch vehicle commit-to-launch decisions include an assessment of the upper-level (UL) 
atmospheric wind environment to determine the vehicle’s predicted controllability and structural integrity 
during ascent. These assessments occur at predefined times during the launch countdown based on the most 
recent measured wind data. However, the pre-launch measured winds may not represent the wind 
environment during the vehicle ascent. Uncertainty in the UL winds over the time period between the 
assessment and launch can be mitigated by engineering analysis of critical vehicle trajectory variables and 
development of statistical models that provide the required level of protection. Trajectory allowances as a 
function of time interval that reduce the maximum allowable values (“knockdowns”) of selected critical 
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variables are made available during the countdown to ensure that the launch decision takes into account the 
uncertainty the wind profiles impose on trajectory assessments. Ideally, these “knockdowns” are derived 
using a large sample of measured wind profiles that represents the ascent wind climatology at the launch 
site. Without these historical data, theoretical statistical models must be used, which can misrepresent 
launch availability. Using an over-conservative model could result in overly restrictive levels of protection 
hence reducing the launch availability. Conversely, using an under-conservative model could result in 
launching into winds that might damage or destroy the vehicle. 

Historical databases were collected from three locations: National Aeronautics and Space 
Administration’s (NASA’s) Kennedy Space Center co-located on the United States Air Force’s (USAF’s) 
Eastern Range (ER) at the Cape Canaveral Air Force Station, USAF’s Western Range (WR) at Vandenberg 
Air Force Base, and NASA’s Wallops Flight Facility (WFF). These historical databases consist of a certain 
number of wind pairs, where two wind profile measurements spaced by the time period of interest define a 
pair. Sets of temporal wind pair databases for five time intervals (0.75, 1.5, 2, 3, and 4 hours) were 
generated from historical data at each location. 

II. Data Sources 

Multiple UL wind measurement instrumentation systems were utilized to maximize sample size in the 
databases for each location. Wind profiles could either be from measurements by a rawinsonde balloon, a 
high-resolution Jimsphere balloon or multiple vertically pointing Doppler Radar Wind Profiler (DRWP) 
systems. The rawinsondes are lofted by latex balloons and transmit the data back to a ground-based 
receiving system. Rawinsondes typically reach 30.5 km (100 kft) before the balloon bursts. Output from 
rawinsonde data is usually presented in pressure levels, which corresponds to unevenly spaced altitude 
levels. To use rawinsonde data for vehicle response assessments, data are linearly interpolated to 30.5 m 
(100 ft) altitude levels to fill in gaps where wind data are not reported. Another balloon-based system, used 
only at the ER and WR for space vehicle support, makes high-resolution 30.5 m (100 ft) wind 
measurements through the use of a specially designed balloon known as a Jimsphere (Wilfong et al. 1997). 
There are two types of high-resolution wind measurement systems that can track Jimsphere ascent. One 
system uses ground-based radar to track a Mylar-coated Jimsphere balloon while another system, known as 



Real Time (RT) Automated Meteorologieal Profiling System (AMPS) High Resolution Flight Element 
(HRFE), uses Global Positioning Satellite teehnology and a clear Jimsphere balloon to track ascent (Divers 
et al. 2000, Adelfang 2003, Wilfong et al. 2000). The AMPS system can also track ascent of a rawinsonde 
balloon, known as a Low Resolution Flight Element (LRFE). The Jimsphere balloon itself is more rigid 
than a rawinsonde, plus it contains roughness elements to reduce self-induced oscillation during ascent 
(Wilfong et al. 1997). The Jimsphere also contains a vent valve in order to maintain a constant volume as 
the balloon ascends. However, maintaining constant volume limits the altitude range the balloon can 
achieve. A Jimsphere can typically reach between 16.7-18.3 km (55-60 kft) (Wilfong et al. 1997). 

Vertically pointing DRWP systems are ground-based instruments that transmit and receive electronic 
pulses that can be converted to wind speed and direction. The DRWP transmitted frequency and antenna 
size dictates the altitude range sampled and the sampling interval. The ER has a 50-MHz and five 915-MHz 
DRWPs that, when their measurements are spliced together, can generate a wind profile from roughly 
0.183-18.3 km (0.60-60 kft) (Barbre 2013). Data from the spliced profile are interpolated to a 15.2 m (50 ft) 
altitude interval (Barbre 2013). Unlike balloon-based systems, the DRWP operates continuously, with the 
50-MHz DRWP reporting measurements approximately every 5-mins and the 915-MHz DRWP reporting 
measurements approximately every 1 5-mins. Both DRWP systems produce wind profiles at vertical 
resolutions acceptable for launch vehicle assessments. These attributes yield orders of magnitude more 
profiles compared to balloon profiles available for developing temporal wind pair databases. The WR also 
has an equivalent set of DRWPs but due to time limitations data from both systems could not be integrated 
to generate vertically complete profiles. 

A. NASA Wallops Flight Facility 

Rawinsondes provide the only source of wind data at WFF, and two databases of rawinsonde profiles 
from WFF were obtained. The first was from the National Climatic Data Center (NCDC) Integrated Global 
Radiosonde Archive (IGRA) (Durre et al. 2006) database for the October 1963 through January 2000 
period of record (POR). The IGRA data for WFF consists of balloons released from the National Weather 
Service (NWS) location near WFF. The other database of rawinsondes was obtained directly from WFF, 
has a POR of February 2000 through January 2013, and consists of rawinsondes released at the NWS site 
and at WFF in support of mission operations. The IGRA database includes the rawinsonde data that was 



directly obtained from WFF personnel, which implies that no reason exists to include the IGRA data post 
December 1999 in the WFF wind pairs database generated. 

B. USAF Western Range 

Archived data from rawinsondes and Jimspheres were available for developing the WR wind pair 
databases. The data came from three sources: IGRA rawinsondes from January 1965 through January 2013, 
WR rawinsonde profiles from February 2008 through April 2012, and a Jimsphere database from January 
1965 through September 2001. The time overlap between the IGRA and WR rawinsonde databases was 
necessary because the IGRA database contains only WR rawinsonde data. Whereas, the WR rawinsonde 
archive contains real-time AMPS (RTAMPS) LRFE wind data in addition to rawinsonde data that is 
archived in the IGRA database. As stated earlier, the WR wind pair databases do not contain DRWP 
measurements because the DRWP systems required extensive time to process. 

C. USAF Eastern Range 

Data from the ER DRWP systems provide the largest sample size and are the sole source used for 
database development. The spliced ER DRWP profde database has a POR of April 2000 through December 
2009. This POR results from the availability of quality controlled (QC’d) data for both the 50-MHz and 
915-MHz DRWP (Barbre 2013). No rawinsonde or Jimsphere data were used because adding these data 
would have only increased the sample size of the ER database by 0.5%. 

Ill, Data Processing and Quality Control Procedures 

Extensive QC of wind profde data were required to remove suspect data in individual profiles as well 
as in profile pairs. Automated and manual QC checks were applied on the data from each measurement 
source. The automated QC checks differed between the measurement sources and consisted of general and 
study-specific checks. The latter checks were necessary because the general QC checks only evaluated 
data in single profiles and did not check consistency within a profile pair. A profile was rejected if it failed 
any given QC check. The following sections present details of the QC checks for each measurement 


system. 



A. Rawinsonde 


Rawinsonde data from all sources went through a two-step QC process. The rawinsonde data obtained 
from IGRA already were processed through a set of QC checks (Durre et al. 2006). An additional set of 
more stringent QC checks were performed on all rawinsonde data based on a manual review of the IGRA 
data. The development process applied the following QC checks to each individual profde: 

At least ten altitudes that contain either wind or thermodynamic data must exist. 

Vector differences between vertically adjacent wind measurements must be less than 

51.4 ms ' (100 kt). 

The mean wind speed over the entire profde must be less than 51.4 m s ' (100 kt). 

Difference magnitudes between vertically adjacent temperature measurements in the lowest 

3.048 km (10 kft) of the profde must be less than 20°C (68°F). 

All heights must increase. 

The minimum altitude must be positive. 

Dew points corresponding to temperatures must be at least -60°C (-76°F). 

The development process first extracted wind pairs containing profiles that passed all initial QC 
checks. However, several more QC checks were applied on the data to check for data quality between 
profile pairs. This process removed a wind pair that contained one or two profiles that failed a QC check. 
The following QC checks were implemented to remove suspect wind pairs: 

Wind data with each profde in the pair must reach a minimum of 6.095 km (20 kft). 

Wind component change between adjacent altitudes (vertical wind shear) must not exceed 

0.15 s-'. 

More than 50% of wind data must exist in both profiles. 

Duplicate pairs were removed. 

Each wind pair was manually inspected for erroneous data. 

The minimum altitude requirement of 6.095 km (20 kft) was based on the minimum altitude required 
in order to perform certain launch vehicle assessments but insufficient for launches to altitudes of 
maximum aerodynamic pressure. Rejecting profiles not containing at least 50% of the possible wind data 



eliminated the potential of having an artifieially large temporal wind ehange at the same altitude over the 
two profiles due to a large data interpolation in one of the profiles. A pereentage eheek as opposed to an 
altitude range eheek was used due to the IGRA data being reported in pressure levels, whieh results in 
unequal altitude intervals and large data gaps inherently existing in valid profiles. This eheek also removed 
profiles that eontained exeessive vertieal wind shear resulting from interpolating data over a large altitude 
interval. As part of the manual QC proeess any profiles that exhibited large interpolations not eonsistent 
with the other profile were removed. The duplieate pair QC eheek was neeessary for the WR wind profiles 
due to using multiple sourees of wind profile databases with overlapping FOR. Several additional 
proeessing and QC eheeks for the WR data were performed due to the overlapping FOR. 

Frofiles from both sourees were merged into a single subset and sorted temporally. 

Frofiles <15 minutes apart were grouped and the profile reaehing the highest altitude was 
ineluded in the database. 

Unique profiles were then merged with the existing IGRA and RTAMFS database. 

The last step entailed manually inspeeting eaeh pair. The development proeess implemented this step 
after a review of the maximum wind eomponent ehange and probability distributions, independent of 
altitude, for eaeh pair time interval. Temporal wind ehange analyses have shown that wind ehange extremes 
are typieally eorrelated to time separation, where longer time intervals imply larger extreme wind ehange 
magnitudes (Johnson, 2000). However, the WFF and WR wind eomponent ehange at probability levels 
greater than 95% in the 0.75 and 1.5-hour pairs were -50% greater than the eorresponding maximum wind 
ehange at the same probability level in the 2, 3 and 4-hour pairs without manual QC. Manual inspeetion of 
the WFF and WR wind pairs for all time intervals revealed these questionable wind ehange values were 
assoeiated with profile pairs oeeurring around 0000 UTC in data obtained from the IGRA database (Durre 
et ah, 2006). This eharaeteristie appeared in the 2, 3 and 4-hour wind pairs for both WFF and WR. The 
differenees observed in the questionable wind profiles seemed more eharaeteristie of diumal-seale wind 
ehange as opposed to short-time period wind ehange, whieh led to questioning the time stamp of eaeh 
profile in the pair. An independent souree provided a eomparison to events of questionable wind ehange. 
Data from the National Centers for Environmental Fredietion (NCEF) North Ameriean Regional 
Reanalysis projeet at the time period of interest were evaluated to determine if a large gradient in the winds 



existed over the time period. The majority of these eases did not eorroborate with the NCEP data and, as a 
result, the wind pairs oeeurring near 0000 UTC that eontained wind profdes from the IGRA souree were 
removed from the generated database. The resultant temporal wind ehange distributions were better 
eorrelated as a funetion of time separation. 

B. Jimsphere 

The WR wind pairs inelude wind profdes from both systems that use the Jimsphere balloon. Jimsphere 
wind profiles were generated during launeh vehiele operations and were manually QC’d by teehnieians 
prior to distribution to launeh vehiele operators. The manual QC eheeks were performed to remove suspeet 
data for use in flight vehiele assessments (Divers et ah, 2000). Additional automated QC eheeks were 
performed on the data prior to eombining with the rawinsonde data to determine where wind pairs existed: 

A profile was removed if its lowest altitude exeeeded 122.2 m (400 ft) or if the profile 
eontained any deereasing altitudes. 

All variables (altitude, wind speed, and wind direetion) were removed if at least one variable 
was missing. 

All data were removed after the first altitude eontaining missing data. 

Any linear wind eomponent interpolations at the top of the profile were removed. 

The seleeted wind pairs for the WR ean be made up of two Jimspheres, two rawinsondes, or a 
Jimsphere and a rawinsonde. The issue with the Jimsphere/rawinsonde eombination is that a differenee 
exists in the smallest resolvable wavelengths between these two wind profiles due to their sampling 
intervals. The small-seale wavelengths were removed through a filtering algorithm in order to maintain an 
equivalent effeetive vertieal resolution between the rawinsonde and Jimsphere systems (Wilfong et ah, 
1997). A 244.3 m (800 ft) filter was applied to the Jimsphere based on a power speetrum analysis of the 
rawinsonde data shown in Fig 1. Filtering the Jimsphere data was neeessary to use wind profiles from 
either system interehangeably in assessing wind affeets on vehiele performanee (Wilfong et al. 1997). The 
RTAMPS mean normalized power speetral density likely eontains additional noise in the 152-610 m (500- 
2000 ft) wavelength range. The additional noise was not removed from the RTAMPS data eontained in the 
database sinee filtering the Jimsphere data to ~610 m (-2000 ft) would remove valid speetral eontent that is 
neeessary to assess wind affeets on vehiele performanee. 



C. DRWP 


The ER wind pairs consist of profiles from the DRWP systems, which are designed to operate 
continuously with limited manual QC processing. The algorithms and methodologies in Barbre (2012) and 
Lambert et al. (2003) detail the QC process on the 50-MHz DRWP (D-50) data and the 915-MHz DRWP 
(D-915) automated QC process, respectively. Additional QC on the D-915 data included removing profiles 
with duplicate timestamps, filling temporal data gaps greater than 15 minutes, and checking for correct 
altitude progression. The resultant POR of temporal overlapping QC’d D-50 and D-915 data extends from 
April 2000 to December 2009. 

The next step after data QC entailed creating a single profile by splicing the two DRWP profiles at the 
altitude where the top of the D-915 profile, 6,100 m (20,013 ft) at most, and the bottom of the D-50 profile, 
approximately 2,500 m (8,202 ft), meet or overlap. Before generating the spliced DRWP profiles, the 
temporal and spatial (vertical) criteria are determined and applied to all DRWP profiles because the 
individual profiles had to match in both domains before splicing (Barbre 2013). The D-50 archive contains 
measurements at 150 m (492 ft) intervals every five minutes prior to an instrument upgrade in 2004 and 
145 m (475 ft) intervals every three minutes thereafter. The D-915 archive contains measurements at 101 m 
(331 ft) intervals roughly every 15 minutes. The spatial criterion for the blended profiles is 50 m (164 ft) 
with the temporal criterion based on the time interval of the D-50 measurements. 

An algorithm was applied to the data to splice the D-50 and D-915 profiles into one using each of the 
five D-915s. First, for each D-50 timestamp, the algorithm found the closest corresponding timestamp in 
each of the D-915 profiles. Next, the algorithm interpolated the D-50 and D-915 profiles to 50 m (164 ft) 
spacing in the altitude range 0.100-18.6 km (328-61,024 ft), which corresponds to the lowest observation of 
the D-915s and the highest observation of the D-50. Then, the algorithm counted the number of data gaps 
from each instrument, flagged 300 m (984 ft) data gaps with the D-50 and 500 m (1640 ft) data gaps with 
the D-915, and linearly interpolated wind components within un- flagged gaps. Both profiles contained data 
placeholders at the same altitudes, with a transition region existing between the two profiles typically at 
altitudes around 2,000-3,000 m (6,562-9,843 ft). The algorithm then spliced the two profiles using a 
methodology that varied slightly based on the data coverage of the two profiles within the transition region. 
If a D-915 profile overlapped the D-50 profile, then the algorithm combined the D-50 and D-915 wind 



components within the transition region using a weighting scheme that provided greater weight to the 
D-915 at lower altitudes and greater weight to the D-50 at higher altitudes. If a D-915 profde did not reach 
the D-50 profile, then the algorithm linearly interpolated the winds between the highest altitude of the D- 
915 profile and the lowest altitude of the D-50 profile, provided the QC algorithm did not flag the gap. 
Each spliced profile contained winds exclusive to the D-915 below the transition region, derived winds 
within the transition region, and winds exclusive to the D-50 above the transition region. Splicing the D-50 
and individual D-915 profiles produced up to five DRWP wind profiles at a given timestamp in the archive, 
depending on the QC process and how many D-915s were operating. 

The individual spliced profiles were then combined to generate a single composite DRWP profile 
representing the wind environment at a given timestamp. The individual spliced profiles only differed 
below the lowest altitude of the D-50 profile. An algorithm was developed that generated composite winds 
at each altitude starting at 150 m (492 ft) and ending at 18,450 m (60,532 ft). The algorithm omitted data at 
the lowest D-915 altitude and the highest D-50 altitude due to ground clutter effects and questionable 
shears, respectively. The algorithm first computed a mean reference wind using the individual spliced 
profiles with valid winds at each altitude. Next, vector differences between the winds from each of the 
individual spliced profiles and the reference wind were used to derive weights corresponding to each of the 
profiles. Summing the product of the weights and wind components from all individual profiles produced 
the composite wind at each altitude. Above the lowest reporting altitude, the algorithm computed the 
reference wind as the mean of the reference wind described above and the composite wind at the previous 
altitude. 

A subsample of the DRWP archive was produced according to specified guidelines. The wind profile 
had to contain data at all altitudes from 250-6,096 m (820-20,000 ft) and the wind components were 
linearly reduced from altitudes below the lowest reporting altitude, which ranged from 150-250 m (492-820 
ft), to no wind at 0 m (0 ft). Profile tops extended as high as 18,450 m (60,532 ft). 

IV. Wind Pair Development and Analysis 

The criteria to select pairs, the resultant number of wind pairs at each location, statistical analyses of 
the sample sizes, and distributions of extreme wind changes are presented in the following sections. 



A. Criteria to select pairs 

Constraining the pair selection to the exact time spacing with the balloon-based WFF and WR profiles 
limits the number of pairs since balloons are released infrequently. Therefore, for each pair the time range 
was expanded by +/- 15 minutes to increase the wind pair sample size. For example, profiles spaced 
between 2.75 to 3.25 hours were treated as 3-hour pairs. For the ER, two profiles defined a pair if the 
desired time separation of the pair +!- two minutes separated the profiles’ timestamps. For example, a 0.75- 
hour (45-minute) pair has two profiles spaced anywhere from 43-47 minutes apart. The pair selection 
process used a two-minute window because a large number of DRWP profiles existed and at least three 
minutes existed between adjacent DRWP profiles. 

B. Sample Size 

Table 1 presents the number of pairs at each time interval and location. Depending on the time interval 
of interest, the WFF, WR, and ER databases contain 54-127 pairs, 366-548 pairs, and 260,878- 
297,490 pairs, respectively. The disparity in the magnitude of samples at the ER is due to the continuous 
operation of the DRWP versus the discrete measurements from the balloon systems used at WR and WFF. 
The WR’s history of supporting space launch operations that require frequent balloon releases attributes to 
the difference between the number of WR and WFF pairs. 

C. Statistical Analysis 

The most frequent application of wind pair databases is to apply the empirical maximum zonal (u) and 
meridional (v) wind change components of each profile into a persistence assessment to determine the 
effects of wind change over a specific time period on vehicle performance (Smith et al. 1992). Therefore, a 
large sample size must exist in order to capture the largest range of maximum wind change possible. 
Several analyses were conducted to determine how well the sample population at each location 
characterized the wind change extremes. 

The analyses results quantify the distribution and the confidence bound (CB) in the empirical 
maximum wind change from the various sample sizes of each pair set. Extreme wind change population 
distributions are usually non-Gaussian (Merceret 1997), so an extreme theoretical probability function was 
used to fit the data. The generalized extreme value (GEV) probability distribution function (PDF) (Coles 



2001, Kotz and Nadarajah 2000) provided a good fit of the extreme u- and v-eomponent wind ehanges in 
eaeh pair up to roughly the 99th percentile level. The GEV PDF is expressed by: 

y = f(x\ k,ii, g) = exp [^1 + fc [^1 + fc ^ (1) 
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where x represents each value in a distribution of wind changes, and k, fx, and a denote the scale, shape, and 
location parameters, respectively, of the GEV estimate. Using the results from the GEV, 95% CB at various 
percentile levels were calculated using the Asymptotic Distribution of Percentiles (ADP) method 
(DasGupta 2008). The ADP equation is a function of the CB, sample size and percentile level of interest. 
The analysis uses the 95% CB as a conservative approach to assess the range of extreme wind change for 
selected percentile levels. 

Representing the data annually in the distribution plots for all sites was necessary due to the limited 
sample populations from WR and WFF. Seasonal and monthly variability does exist in wind component 
distributions (Johnson 2000). When using the data for engineering assessments, ideally, data should be 
grouped by season or month to ensure the “knockdown” value is based on a homogenous sample 
population. However, a database subset representing a given month or season must contain a sample size 
large enough to characterize the environment to the desired probability level. The ER database, which 
contains enough samples for monthly assessments, was compared to the 2013 Cape Canaveral Air Force 
Station Range Reference Atmosphere (RRA) (Bums 2013). The RRA contains statistics of wind variables 
derived from rawinsonde balloon measurements at the ER. Statistics are reported as a function of altitude 
for each month and annually. The comparison first computed the statistics listed in the RRA from the ER 
database, and then overlaid results from both sources. The annual results compared very well with some 
minor differences for an individual month as shown in Fig. 2. These differences are likely caused by the 
different sampling characteristics between the DRWP systems and the balloon-based systems that were 
used to derive the RRA. 

Annual distribution plots of the maximum change in wind component magnitude are presented in 
Figures 3-6 to illustrate differences in the CB as a function of sample population. For each plot, the 
cumulative probability, drawn from the probability density function (Wilks 2006), is along the y-axis and 



the magnitude of the wind component’s change is along the x-axis. The sample size of the pairs is 
correlated to the width of uncertainty at the 95% CB for the highest percentile levels in the sample 
population. As the sample size decreases the width of uncertainty at the 95% CB increases. Figures 3 and 4 
show the variation in the width of uncertainty at the 95% CB. In addition, a small probability density at a 
selected percentile level increases the width of uncertainty. Figure 3, the WFF 2-hour pairs plot, shows this 
attribute - where the 95% CB in the v-component change significantly exceeds the CB in the 
corresponding u-component change even though the sample sizes for both u- and v-changes are the same. 

The width of uncertainty in the CB for all the ER pairs is small due to the large sample size. An 
example of an ER pair fit to the CB is shown in Fig. 4. The deviation of the CB from the empirical 
distribution above the 95th percentile level is an artifact of the CB being calculated from the GEV function, 
which does not fit the empirical distribution well. Flowever, the poor fit is not an issue since the sample size 
is large enough to justify using the empirical percentiles for almost any flight vehicle assessment. 

For the WFF and WR samples, respectively, the 95% CB range of uncertainty increases as the sample 
size decreases. The WR 95% CB range of uncertainty at the sample size’s maximum empirical probability 
level was approximately 15.4 m s ' (30 kt) for both wind components in all the pairs except for 4 hours. The 
range of uncertainty is -41.2 m s ' (80 kt) as shown in Fig. 5. Because of the large uncertainty at the 
extreme empirical percentile in the 4-hour pairs, another approach was applied to quantify the confidence 
of the empirical wind change data. This approach uses a function from Smith and Adelfang (1998) that 
approximates the probability level of a sample population with a specified sample size to a probability level 
of the universal population. The function makes no assumption to the form of the probability distribution 
function of the wind change and is defined as: 

= 1 + [(n - 1) - P” (2) 

where P„ is the probability that the sample contains the universal population at the sample probability Ps 
and the sample size, n. Stated another way; a certain sample size is required to be Pu percent confident the 
sample contains the Ps value of the universal population. Table 2 presents the confidence level of the 
universal population for various sample probability levels based on the sample size in each WR wind pair 
interval. For the 366 4-hour wind pairs, there is 88.1% confidence that the pairs contain the 99th percentile 
of wind change during this time interval. The confidence level exceeds 90% for the other four time periods. 



These samples are large enough for most vehiele performanee applieations; however, a low eonfidenee 
exists that these samples eapture wind ehanges at extreme (e.g., > 99th pereentile) levels. 

The WFF samples eontain the smallest number of pairs of the three loeations. Due to the small sample 
sizes for eaeh time period, the 95% CB for the observed wind ehange extremes (~>20.6 m s ' (40 kt)) have 
a large uneertainty, whieh is more pronouneed for the v-eomponent (Fig. 5). At eaeh time period the 95% 
CB for the v-eomponent wind ehange range is at least 20.6 m s ' (40 kt). The maximum 4-hour 
v-eomponent wind ehange of 37.0 m s ' (74 kt) exists at the 98th pereentile level in the sample population’s 
distribution. The 95% CB at the 98th pereentile level ranges from 20.7 to 45.9 m s ' (40.2-89.3 kt) as shown 
in Fig. 6. 

Table 3 presents eonfidenee levels of the universal population for various sample probabilities based 
on the WFF sample size. A 16.9% eonfidenee exists that the 4-hour pairs eontain the 99th pereentile of all 
wind ehanges during this period. The eonfidenee levels range from 10-36% at the 99th pereentile for the 
other pairs. Due to the low eonfidenee that the sample eontains extreme wind ehanges and large uneertainty 
in the 95th eonfidenee intervals at probability levels above 95%, the reeommendation is to apply the 
extreme 4-hour wind eomponent ehange for all time ehange intervals of interest in vehiele performanee 
evaluations. Applying this reeommendation produees more eonservative results for shorter time periods, 
while inereasing the likelihood of under-eonservative results as the time period approaehes 4-hours. 

Table 4 presents eonfidenee levels of the universal population for various sample probabilities based 
on the ER sample size. The eonfidenee level is 100% for all time periods. 

V, Conclusion 

Temporal UL wind pair databases were generated for ineorporation into eommit-to-launeh deeisions 
based on UL wind assessments. Databases for five time intervals (0.75, 1.5, 2, 3, and 4 hours) at the USAF 
ER and WR, as well as NASA’s WFF were generated through use of historieal data at eaeh loeation. 
Multiple sourees that measure UL atmospherie winds at the requested sites were used for this task. 
Databases were eompiled using wind profiles from rawinsonde, Jimsphere, and DRWP systems. Extensive 
QC eheeks were applied on the data to remove unaeeeptable profiles. Statistieal analyses of the resultant 
wind pairs from eaeh site were performed to determine if the observed extreme wind ehanges in the sample 



pairs represent extreme temporal wind ehange. The resultant ER wind pair databases yielded sample sizes 
that eharaeterize the extreme wind ehange environment and no restrietions on the usage exist. The WR 
wind pair database sample size is large enough for vehiele performanee assessments up to the 99th 
pereentile level. However, due to the small sample size for eaeh wind pair time period at WFF, low 
eonfidenee exists that the observed extremes in eaeh time period eharaeterizes the extreme wind ehange 
environment. Therefore, for any vehiele performanee applieations at WFF, the reeommendation is to apply 
the extreme 4-hour wind ehange values for all time ehange intervals of interest. 

Future work would inelude inereasing the number of WR wind pairs by adding data from the WR 
DRWP systems into the WR temporal wind pair databases. This process would include, at the minimum, 
QC of the individual 50-MHz and 915-MHz DRWP wind profiles. Acceptable wind profiles from each 
DRWP system would be spliced to generate vertically complete wind profiles and available pairs would 
then be incorporated into the existing databases. 
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Table 1 Sample size of wind pairs at each location. 


Time Interval 
(hours) 

ER 

WR 

WFF 

0.75 

273,265 

435 

78 

1.5 

260,878 

401 

54 

2 

297,490 

548 

75 

3 

273,189 

508 

127 

4 

276,108 

366 

74 

TOTAL 

1,380,930 

2258 

408 


Table 2 Confidence levels of the nniversal popnlation for arbitrarily selected sample probability 
levels and the WR sample size for each wind pair time interval (Smith and Adelfang 1998). 


Sample 

Probability 

Time Interval (Sample Size) 

0.75 hours 

1 .5 hours 

2 hours 

3 hours 

4 hours 

(435) 

(401) 

(548) 

(508) 

(366) 

0.500 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.750 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.900 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.950 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.990 

0.9318 

0.9102 

0.9734 

0.9628 

0.8813 

0.995 

0.6400 

0.5960 

0.7592 

0.7215 

0.5466 

0.999 

0.0711 

0.0617 

0.1050 

0.0925 

0.0525 















































































Table 3 Confidence levels of the nniversal popnlation for arbitrarily selected sample probability 
levels and the WFF sample size for each wind pair time interval (Smith and Adelfang 1998). 


Sample 

Probability 

Time Interval (Sample Size) 

0.75 hours 

1.5 hours 

2 hours 

3 hours 

4 hours 

(78) 

(54) 

(75) 

(127) 

(74) 

0.500 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.750 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.900 

0.9973 

0.9763 

0.9965 

0.9999 

0.9962 

0.950 

0.9065 

0.7592 

0.8944 

0.9886 

0.8900 

0.990 

0.1836 

0.1018 

0.1729 

0.3629 

0.1693 

0.995 

0.0584 

0.0301 

0.0545 

0.1332 

0.0532 

0.999 

0.0028 

0.0013 

0.0026 

0.0073 

0.0025 


Table 4 Confidence levels of the nniversal popnlation for arbitrarily selected sample probability 
levels and the ER sample size for each wind pair time interval (Smith and Adelfang 1998). 


Sample 

Probability 

Time Interval (Sample Size) 

0.75 hours 

1.5 hours 

2 hours 

3 hours 

4 hours 

(273,265) 

(260,878) 

(297,490) 

(273,189) 

(276,108) 

0.500 

1 

1 

1 

1 

1 

0.750 

1 

1 

1 

1 

1 

0.900 

1 

1 

1 

1 

1 

0.950 

1 

1 

1 

1 

1 

0.990 

1 

1 

1 

1 

1 

0.995 

1 

1 

1 

1 

1 

0.999 

1 

1 

1 

1 

1 
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Fig 2. Comparison of ER DRWP profiles contained in 4-hr pairs database with the 2013 CCAFS 
RRA database statistics for February (a) monthly mean u- and v-component (solid line - CCAFS 
RRA; dashed line ER DRWP), June (b) monthly mean u- and v-component (solid line - CCAFS 
RRA; dashed line ER DRWP) and Annual (c) mean u- and v-component (solid line - CCAFS RRA; 
dashed line ER DRWP). 
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with Generalized Extreme Value Fit and 
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b) 

Fig 3. WR maximum wind change from the 4-hour wind pairs with 95% CB for the u-(a) and v- 
eomponent (h) wind changes. The magnitude of the wind component change is on the x-axis and 
prohahiUty is on the y-axis. The number (n) of pairs in the analysis is 366. 


ER Maximum 4-hr u-component Change Magnitudes 
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ER Maximum 4-hr v-component Change Magnitudes 
with Generalized Extreme Vaiue Fit and 
95% Confidence Bounds, n =276108 



b) 

Fig 4. ER maximum wind change from the 4-hour wind pairs with 95% CB for the u-(a) and v- 
component (h) wind changes. The magnitude of the wind component change is on the x-axis and 
prohahility is on the y-axis. The number (n) of pairs in the analysis is 276,108. 
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b) 


Fig 5. WFF maximum wind change from the 2-hour wind pairs with 95% CB for the u-(a) and v- 
component (h) wind changes. The magnitude of the wind component change is on the x-axis and 
prohahiUty is on the y-axis. The number (n) of pairs in the analysis is 75. 
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WFF Maximum 4-hr v-component Change Magnitudes 
with Generalized Extreme Value Fit and 
95% Confidence Bounds, n = 74 



Fig 6. WFF maximum wind change from the 4-hour wind pairs with 95% CB for the u-(top) and 
component (h) wind changes. The magnitude of the wind component change is on the x-axis and 
prohahiUty is on the y-axis. The number (n) of pairs in the analysis is 74. 



